Unstable Hadrons in Hot Hadron Gas: in Laboratory, and in the Early Universe 
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We study kinetic master equations for chemical reactions involving the formation and the natural 
decay of unstable particles in a thermal bath. We consider the decay channel of one into two 
particles, and the inverse process, fusion of two thermal particles into one. We present the master 
equations for the evolution of the density of the unstable particles in the early Universe. We obtain 
the thermal invariant reaction rate using as an input the free space (vacuum) decay time and show 
the medium quantum effects on n + n p reaction relaxation time. As another laboratory example 
we describe the K + K <j> process in thermal hadronic gas in heavy ions collisions. A particularly 
interesting application of our formalism is the ty° o 7 + 7 process in the early Universe. We also 
explore the physics of ir^ and fj, freeze-out in the Universe. 

PACS numbers: 95.30.Cq, 52.27.Ny, 24.10.Pa 

I. OVERVIEW 

A. Particles in the Universe for T > 5 MeV 

This study began with the question: At what temper- 
ature in the expanding early Universe does the reaction 



7T <-» 7 + 7 



(1) 



'freeze' out, that is the tt° decay overwhelms the produc- 
tion rate and the yield falls away from chemical equilib- 
rium yield. Because the tt° life span (8.4 10 -17 s) is rather 
short, one is tempted to presume that the decay process 
(arrow to the right) dominates. However, there must 
be a detailed balance in the thermal bath: the produc- 
tion process (arrow to the left) in a suitable environment 
must be able to form ir° with strength corresponding to 
the decay process lifespan. 

We demonstrate here that the 7T° production and equi- 
libration relaxation time is of the same order of magni- 
tude as the lifespan of ir° in the post-quark-gluon-plasma 
hadronization Universe, T < 200 MeV. The point is that 
the 7T° life span is much shorter than the Universe expan- 
sion time (inverse expansion rate) 1 / H Q : 
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where g* is the number of degrees of freedom. m p i — 
1.2211 10 19 GeV is the Plank mass. Figure [1] compares 
the 7T° production-equilibration time (blue solid line) 
with the Universe expansion time 1/H dashed (green) 
line]. We see that tt° equilibration time is much shorter, 
by 14 orders of magnitude at T = 10 MeV, compared to 
the Universe expansion time constant. 

The reason for this is that in thermal equilibrium the 
photon density remains high [dash-dotted (red) line in 
figure [5] also for relatively small T. Thus there is a 
small non-negligible probability of finding high energy 
photons capable to produce ir°, whose density at low T 
is very small (solid blue line in figure [2]). The ir° pro- 
duction has enough time to equilibrate with the decay 



FIG. 1: (Color online) it equilibration time [solid (blue) line] 
and Universe expansion time 1/H as functions of temperature 
[dashed (green) line]. 



process. Therefore the 7r° density does not freeze out, 
but decreases with decreasing ambient temperature of 
the expanding Universe, all the time remaining in chem- 
ical equilibrium with the photon abundance. 

Let us recall how the Bose distribution describes n°- 
density: 



eg 
n Jj = 



d 3 p 



1 



(2tt) 3 e"-P/ feT - 1 ' 



(3) 



Here it M is the four- velocity of the observer with reference 
to the heat bath rest frame, p M is momentum four- vector 



!>' = [—,P 



E = y/ 'p 2 + m 2 



(4) 
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FIG. 2: (Color online) Thermal equilibrium density as a func- 
tions of temperature for: 7 [dash-dot (red) line], 7r° [solid 
(blue) line], ^ pair [dashed (green) line], and nucleons p + n 
[dotted (black) line 3]. 

of the particle considered: a similar expression applies 
for photons, which have two fold spin degeneracy, and 
Tn — y 0. The resulting 7r° density falls exponentially when 
kT < m^oc 2 (henceforth units are chosen such that h — 
c = k = 1). However this density remains high compared 
to the nucleon's density in the Universe (dotted (black) 
line in Fig. [2j taken from Q), down to a temperature 
of about 6 MeV. This is the lower T-limit of validity 
in our present study, as we consider particle production 
reactions in a particle- antiparticle symmetric Universe. 

Some of the results we derive here were presented in Q 
without a derivation: there we considered a laboratory 
e + e _ 7 plasma and postponed the theoretical and ana- 
lytical details. Here we evaluate the reaction relaxation 
time for reactions involving two particles fusing into one 
particle, and/or particle decay into two, Eq. ([T]), and re- 
late this to the lifespan of decaying particle in vacuum. 
To complement this in chapter HVl subsection (c) we also 
consider ir^, which can be equilibrated by the reaction: 

7T° + 7T° <-» 7T+ + 7T~. (5) 

and also by reactions involving muons 

7T ± -H- /U ± + ^(^ (tl ), (6) 

We also consider how fast muons are produced in the 
reactions: 

7 + 7 <-» /!+ + jiT , e + + e~ «-> fi + + pT ; (7) 

and show that these particles also do not freeze out down 
to a T value of a few MeV p|. The muon density is 



slightly higher than that pions because of their smaller 
mass; see the dashed (green) line in Fig. [2] Another 
reaction that may influence muon chemical equilibration 
is the decay of one particle to three particles, and the 
reverse reaction, including neutrinos: 

^ <-> e ± + v e {v e ) + v^Vn). (8) 

However in this case the exact influence of medium ef- 
fects on the reaction rate is more complicated and we do 
not consider this reaction in complete detail here. We do, 
however, compare the relaxation times of particle produc- 
tion in all mentioned reactions with the Universe expan- 
sion rate to see if the particle densities stay in chemical 
equilibrium. 

B. Degrees of Freedom in the Universe 

The lifespans of all unstable hadrons and leptons, ex- 
cept for neutrons n, are much shorter than the Universe 
expansion rate for 5 < T < 200 MeV. Here we show 
that, as a result, all unstable particles stay in chemical 
equilibrium, including neutrons which are effectively sta- 
ble on the time scale of expansion. The importance of 
this remark is that we can evaluate the active effective 
degeneracy (degrees of freedom) in the Universe in the 
entire temperature domain including all unstable hadron 
states. In the hadron phase we define the effective de- 
generacy using as reference the Stephan-Boltzmann law, 

9e(T) = ^, *=^, (9) 
where e is the energy density: 

e = f Y,9 i E l f l (p)d 3 p, E t = yfiri[+iP, (10) 
"* i 

whith the sum over all particles present. 

In Fig. [3] we show the degeneracy Eq. ^ as a func- 
tion of T. The dashed (red) line accounts for the pho- 
ton, three families of neutrino and antineutrino, electron, 
positron, muon, and antimuon contributions. The dot- 
dashed (green) line adds pions; the solid (blue) line, all 
hadrons. Pions begin to contribute noticeably to degen- 
eracy at T > 30 MeV. Among hadrons we included all 
light and strange mesons and baryons, up to a mass of 
about 1700 MeV. The finite density of p and n is also in- 
cluded. As noted earlier, this is a more complicated case; 
fortunately the finite baryon density contributes at just 
at a few percent to gE near the hadronization tempera- 
ture, where the particle-antiparticle symmetry is good to 
10 orders of magnitude. 

The boundary between quark-gluon and hadron phase 
(vertical line) is near T/, = 170 MeV. In figure [3] we 
also show degeneracy in QGP for T > 160 MeV (up- 
per lines). The results shown are based on our earlier 
detailed study of QGP properties @- Here, the QGP de- 
generacy is shown for the extreme cases of either strange 
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FIG. 3: (Color online) Effective degeneracy gs in the Universe 
based on the energy density of hadrons, and for QGP, as a 
function of T. See text for more details. 



quark m s = 90 MeV [dashed (purple) line] or 160 MeV 
[solid (turquoise) line] . Since the expansion of the Uni- 
verse is relatively slow compared to expansion of QGP in 
laboratory, heavy and strange quarks also have enough 
time to reach chemical equilibrium density in the QGP 
temperature range presented in graph. 

We see in figure [3] that the effective degeneracy of 
hadrons, while rising fast, is still smaller than the degen- 
eracy in QGP in the domain of phase transformation tem- 
perature near 160-170 MeV. Many heavy hadron states 
may be missing from the experimental tables. Even 
though their individual contribution to the degeneracy 
is decreasing, their number is expected to grow rapidly, 
in accordance with the Hagedorn hypothesis, in which 
hadron mass spectrum diverges exponentially near the 
hadronization temperature. This theoretical exponen- 
tially growing component leads to a smoother transition 
between hadronic gas and QGP, as is qualitatively indi- 
cated in Fig J3] [dotted (black) line]. 



C. Production and decay of unstable particles 

We show here for the first time the detailed derivation 
of relaxation time for reactions involving one-to-two par- 
ticles in the thermal medium, which we considered in Q 
- In the rest frame of the decaying particle to 3 , the 
reaction 



requires that m\ + to 2 < m^, which allows the sponta- 
neous decay process. This is easily seen considering 



(pi +P2) 2 

(mi + m 2 ) 2 + 2(EiE 2 - mim 2 - p\ ■ p%) 



> (mi + m 2 y 



(12) 



In the last inequality we used E\E\ > (mim 2 + Pi -pi) 
which can be reorganized to read {m\P2 — m-ipx) 2 > 
Pi • Pi — P1P2 ■ This is always true since the right hand 
side is always negative, or zero if both vectors are parallel. 
The equality sign corresponds to the case mi + m 2 = 1113, 
where the reaction rate vanishes by virtue of vanishing 
phase space. This text-book exercise shows that the re- 
action Eq. ([T]) is possible when condition Eq. (|12p is sat- 
isfied. 

The constraint Eq. (fT2|l forbids many reactions. For 
example, the hydrogen formation p + e — >H is forbidden, 
as for a bound state, m# < m p + m e . Thus there must 
be a second particle in the final state. The electron cap- 
ture involves either a radiative emission, p + e — >-H+7, or 
a surface/third atom, which picks the recoil momentum. 
The situation would be different if there were "resonant" ' 
intermediate states of relatively long life span with en- 
ergies above the ionization threshold. Such "doorway" 
resonances are available in many important physical pro- 
cesses. 

It is natural to evaluate the rates of the processes of 
interest, Eq (Eq. (fTTj) ) in the rest frame of particle '3', 
boosting, as appropriate, from or to laboratory frame. To 
do this effectively we need the master population equa- 
tions in an explicitly covariant fashion, which is discussed 
in Sec: [TT1 see Ref. [6j. The kinetic equation for time evo- 
lution of number N of decaying particles 3 can be written 



1 dN 3 
V~dT 



T 1 T 2 
T 3 



- 1 



dW- 



3^12 



dVdt 



(13) 



where dW3^i2/dVdt is the decay rate of particle 3 and 
Ti is the fugacity of particle i. Here the number density 
rii of particle i in thermal (kinetic), but not necessarily 
in chemical, equilibrium is given by 



Ni 

fb/fi^UPi) 



(2tt) s 



d 3 Pif b /f(p l ), (14) 
(15) 



n C7 1 e( u -P*-* t *)/T =p r 



/ is the covariant form of the usual Bose or Fermi dis- 
tribution function defined in the rest frame of the ther- 
mal bath, and describes the corresponding quantity in a 
general reference frame where the thermal bath has the 
relative velocity defined by u M . In the rest frame of the 
thermal bath frame we have: 



1,0). (16) 
Pi is the four- momentum- vector of particle i: 



A!+B 2 « C" 3 



(11) 



P 



1" _ 



(Ei,Pi) ■ 



(17) 
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/ij is the chemical potential, which shows the asymme- 
try in particle and antiparticle densities pi = — fy. For 
reactions considered here we have pi ~ 0. This was as- 
sumed in Eqll3l Note that the distribution function / is 
a Lorentz scalar but the spatial density rii is not. 

Particle C 3 attains the chemical equilibrium when the 
following condition among fugacities is satisfied: 

TiT 2 = T 3 . (18) 

This, as expected, is equivalent to the Gibbs condition 
for the chemical equilibrium. In Sec. IIIII we evaluate the 
invariant rate using the vacuum decay time established 
in the rest frame of the decaying particle, and we discuss 
the behavior of the average decay rate of an unstable par- 
ticle in the presence of the thermal bath. In Sec. IIV1 we 
apply our formalism to two examples: 

a) We study the formation and decay rate of the p 
meson through n + tt <-> p in a baryon-free hot hadronic 
gas, where mesons are considered in thermal and chemi- 
cal equilibrium. 

b) We consider the reaction 7 + 7 <H- ir° in the early 
Universe and find that the expansion of the Universe is 
slow compared to pion equilibration, which somewhat 
surprisingly (for us) implies that 7r° is at all times in 
chemical equilibrium (but at sufficiently low tempera- 
tures e.g. 3-4 MeV, the local density of ir° is too low 
to apply the methods of statistical physics) . 

c) We consider the reaction Eq.© as an example of 
the decay of w to fermions, and the reverse reaction, and 
show that 71"^ and p^ are also in chemical equilibrium un- 
til their equilibrium density vanishes at low temperatures 
(about 3-4 MeV) because of the large mass. Also, we dis- 
cuss neutrinos equilibration by way of this reaction. 

d) We study <fi mesons evolution considering the reac- 
tion K + K •<-» <j) in thermal hadronic gas in heavy-ion 
collisions. 

To conclude this overview we draw attention to the 
fact that unlike the one-to-two reaction the two-to-two 
reactions 

Ai + B 2 O C 3 + Di (19) 

have been extensively studied in the past, in the context 
of astrophysics and cosmology 0, Q and heavy ion reac- 
tions 5]. However, the simpler one-to-two situation has 
escaped attention so far, and the adaptation of kinetic 
methods is in detail not trivial given the novel quantum 
and relativistic effects involving particle decay. Aside 
from cosmology implications, the insights gained in this 
study are clearly of relevance to the general understand- 
ing of QGP and hadron gas evolution in relativistic heavy 
ion collisions. For example our present work allows to 
consider the chemical yields arising in reactions such as 
p o 7T7T, 7T° 077, Ao Nir and K + K o cj> [|, Q. 



II. KINETIC EQUATIONS FOR DECAYING 
PARTICLES 

A. Decaying particle density evolution equation 

Consider an unstable particle, say C3, which decays 
into other two particles, 

C 3 -> Ai + B 2 (20) 

in the vacuum. In a dense and high-temperature thermal 
ambient phase, particles Ai and B 2 are present, and the 
inverse reaction, 

A x + B 2 -> C 3 (21) 

can occur, producing the particle we called C 3 . For 
now we assume that the abundance of particle C 3 
changes solely by decay Eq. (|20|) and (thermal) produc- 
tion Eq. (|2ip reactions. The time variation of the number 
of particles N 3 than is controlled by the master equation 

1 dN 3 = dW 12 ^ 3 dW 3 ^ 12 

V dt dVdt dVdt ' 1 ' 

where dWi 2 ^, 3 /dVdt is the production rate per unit vol- 
ume of particle type C 3 via Eq. (|2"Tj) and dW 3 ^i 2 /dVdt 
is the decay rate of particle type C 3 per unit volume. 

A very similar master equation controls the abundance 
of particle A\ and B 2 

1 dN h2 _ dW 3 ^ 12 dW 12 ^ 3 t 

V dt ~ dVdt dVdt [Z6 > 

where the rate i? thor is caused by other reactions influ- 
encing the abundance of particles of type Ai and B 2 . 

As an example, consider the reaction p -h- tttt in dense 
hot matter formed in heavy-ion collisions. Pions can 
be easily created by inelastic collisions of other hadrons 
and thus we have to deal with a multicomponent sys- 
tem involving i? ther when looking at ir abundance, but 
to evaluate p abundance the dominant terms are as in 
Eq.(Eq. (|2"2"|) h We often can assume that i? ther domi- 
nates the yield gains and losses and thus we can use the 
thermal distribution for the particles A\ and B 2l which 
in the example above are pions. 

In the following, we thus assume that particles A\ and 
B 2 are in thermal equilibrium and further, we assume 
that the system is spatially homogeneous. In a thermal 
equilibrium, the dynamical information can be obtained 
from the single particle distribution function / (p) for 
each particle, see Eq. (fl"5|) . / is controlled by two pa- 
rameters, the temperature T and the fugacity T . In this 
paper, we assume that the fugacity T changes over time 
by way of chemical reactions much more rapidly than 
does the temperature T of the ambient thermal bath, 
and thus we can consider reactions at a given constant 
T. This assumption is certainly valid in the domain of 
temperatures we consider, and may fail only at the very 
highest primordial T in the early Universe. 



■5 



B. Decay and production rates 

The thermal production rate dWn^s/dVdt and the 
decay rate of particle 3 under the thermal background 
dW 3 ^\ 2 /dVdt can then be expressed using these distri- 
bution functions for each of the particles involved in the 
reaction. According to the boson or fermion nature of 
the particle Ax, we have to consider different cases. If 
particle A\ is a boson, then there are two cases of the 
decay and production mode, and if the particle C 3 is a 



fermion, it can only decay into a boson and a fermion: 

boson3 i — > bosoni + boson2, (24) 
boson3 4 — > fermioni + fermion . (25) 
fcrmion3 < — > bosoni + fermion . (26) 

Accordingly, the Lorentz invariant transition probabil- 
ity per unit time and unit volume corresponding to the 
process Eq. pTj) is 



dW 12 ^ : 
dVdt 



,9i 



d 3 Pi , ,„ , 92 
/b,/(Ti,pi)- 



d 3 P2 
2E 2 



/6,/(T 2 ,p 2 ) 



d 3 p 3 



1+J(2tt) 3 J 2E 1 J "^ jy l ^ ±j (2tt) 3 J 2E2"" y J 2E 3 {2k)° 

x (2tt) 4 5 4 (pi + P2 - p 3 ) — ]T \(PiP2 \M\ P3 )| 2 (1 ± fb,f(r 3 ,p s )) , 



5152 



(27) 



spin 



where / = 1 for the case of a reaction between two in- 
distinguishable particles A\ and A 2 , and I — if A\ and 
B 2 are distinguishable. The factor 1/(3152) and the sum- 
mation are caused by averaging over all initial (iso)spin 
states. The last factor in Eq. (|2T|) accounts for the en- 
hancement or hindrance of the final-state phase owing 
to the quantum statistical effect, as introduced first by 
Uchling and Uhlenbeck [9(. The upper sign '+' is for the 
case when particle C3 is a boson; the lower sign '— ' when 
it is a fermion. Eq. (|2T[) is manifestly Lorentz invariant 
and therefore it can be used in any frame of reference. 

This rate is related by a detailed balance relation 
to the particle C3 decay rate 



dWu 



dW 3 



dVdt 



dVdt 



-T1T9 



(28) 



Therefore chemical equilibrium T1T2 = T3 corresponds 
to the condition of equal decay and production rates as 
we expected. Using Eq. (f2"5]l . Eq. (f2"2"j) can be written in 
the form of Eq.fEq. (|T5|)). 

Equation (fT3|) can be further simplified by defining the 
decay time in matter of particle i: 



R 



where rate: 



R = 



1 dW : 



3^12 



1 dW, 



12^3 



T 3 dVdt TiT 2 dVdt 



(29) 



(30) 



We show in the next section that this definition has the 
right vacuum limit and that the dynamical equations as- 
sume a particularly simple form. However, the reader 
should observe that other definitions could be considered. 

It is convenient to introduce kinematic reaction times 
in analogy to the dynamic expression Eq. (I2U1) . Doing 



this we cast Eq. (|13[) into the form of an equation for T 3 
* 3 = — T 3 + — T 3 + -(TiT 2 -T 3 ), (31) 

T~T T~S T 3 

where we defined the kinematic relaxation times related 
to the evolution of temperature and entropy 



— = -T 3 g 



1 

t t 
1 

TS 



:i ^ d(n v /(T 3 g*T 3 ))/dT 
dn^/dTg, 



T, 



n n /T 3 d\n(g*VT 3 ) 



dn 7r /dT 3 



dT 



T. 



(32) 
(33) 



We introduced the minus sign above in order to have tt, 
ts > 0. Compared to our earlier definition [6J, we now 
included g* in tt and ts- 

While in principle the values of tt and t$ are unre- 
lated, for a given kinematic stage of system evolution the 
temperature change can be related to entropy change. 
For example, for the radiation-dominated epoch of the 
Universe we have 



T 
T 



R 
R' 



(34) 



In the radiation-dominated Universe the entropy conser- 
vation further implies that 



1 

TS 



0. 



Freeze-out from chemical equilibrium arises for 



tt ~ r 3 



(35) 



(36) 



When tt is smaller than t 3 , that is the kinematic term in 
Eq. pip is important, T3 begins to increase often rapidly, 
as the number of particles '3' is preserved but their den- 
sity decreases owing to dilution in expansion. Because 
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T3 > 1 the multiplicity of particle '3' in the slow de- 
cay is dominant, especially so when the particles 1 and 
2 yields remain in chemical equilibrium by the action of 
other processes. These considerations can be of great im- 
portance in the study of dark-matter particle abundance 
where the lifespan against decay and/or annihilation is 
comparable to the life span of the Universe. We postpone 
further discussion to a more appropriate opportunity. 



III. EVALUATION OF THE INVARIANT 
DECAY (PRODUCTION) RATE 

A. General case 



The vacuum decay width of particle C 3 in its own rest 
frame can be found in textbooks. In our notation, 



J 



2m 3ff3 4 (7 + 1) (2?r) 



ITTirf / ik 5tEl + E2 ~ m ' ] S m ~ ml m3) 



11 1 

r ~~ 2m 3 1 + I J 2E 1 (2nf J 2E 2 (2tt 
I 1 r ^3 



8m 2 3 g 3 (I + 



J— ^2m-p\M\m 3 ) 



spin 



(37) 



r 



Here p = Pi = P2 and = ^/P z + fn% 2 are, respec- 
tively, the magnitudes of the momentum and the ener- 
gies, of the two particles A\ and 5 2 in the rest frame of 
particle C 3 



P 



m 3 ± (mf — m|) 



2m 3 



(to? 



4mg 



(38) 



The magnitude of three-momentum \p\ is of course the 
same for particle A\ and B2 in the rest frame of decaying 
particle C 3 . 

We denote by t 3 the decay rate of particle C 3 in the 
rest frame of the thermal bath in which it is emerged, E 3 



and p 3 are the corresponding energy and the momentum 
of particle C 3 which changes with thermal velocity distri- 
bution. The thermal decay reaction rate per unit volume 
dWz^i+2 1 'dV dt is then obtained by weighting l/r 3 with 
the probability to find the particle at a given momentum 
and introducing the Lorentz factor 7, so that E 3 T 3 /m 3 is 
the decay time of particle 3 with moment p 3 



dW : 



3->l+2 



dVdt 



93 
(27T) 



^3 / d 3 p 3 f bJ (T 3l p 3/ 

71") J 



m 3 1 



(39) 



Comparing Eq. (j3"9")l with Eq. (j2"7) , we conclude that in 
medium, at finite temperature T, the decay rate r 3 of 
particle C 3 in the rest frame of the heat bath is given by 



J 



1 1 



d 3 pi 



t' z 2m 3 1 + I J 2Ex (2tt) 3 J 2E 2 (2tt 



d 2 P2 ,„ >4 ; 4/ , \ 
3 / n rT7r^3 ( 2lT ) ° (Pl+P2~P3) 



x - ^|(piP 2 |Mb 3 )| 2 / & ,/(Ti,Pi)A,/(T 2 ,p 2 )Tr 1 T 2 - 1 exp( w .p 3 /T), 

#3 

spin 



(40) 



which is a Lorentz invariant form, but u ■ p 3 — > E 3 , the 
energy of the particle 3 in the rest frame of the thermal 
bath. 

Using the in-vacuum particle C 3 rest-frame decay time, 
Ea.([37]l. we find that Eq.gDJ takes the form: 



1 1 e E ^ T 

To T 2 



(41) 



The function $(p 3 ) is: 

$ (pa) : 
with 



1 T^e^ 



-60 ± i r -i e (a 2 +bO ± 1 



(42) 



EiE 3 E 2 E 3 pp 3 

«i = ^ a 2 = =r, 6 = = and 

m 3 i m 3 J m 3 J 

£ = cos 9 = cos(p*2 A pi). (43) 
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With this the particle C 3 decay rate per unit volume in 
a thermally equilibrated system is given by 



dVdt 



e 3 /t 



(2tt 2 ) to Jo E 3 T3 W t ± 1 



$(ps), 
(44) 

We were able to evaluate the integral $(^3) analyti- 
cally in the absence of particle- antiparticle asymmetry 
(absence of chemical potentials) , 



IV. EXAMPLES 
A. Hadrons in heavy-ion collisions 

1. Production of p mesons via the p tttv process 

First, we consider an example of p-meson thermal de- 
cay and production in a thermal and chemically equili- 
brated pion bath: 



1 



5 ( e a 1+ a 2 ±T 1 T 2 ) 



In 



(T 2 e- a2 ± e b ) (e ai ± Ti 
(T 2 e- Q 2 ± e ~ b ) (e°i ± Tie 6 



P° O 7T + + 7T , 
P ± ■(^■W ± + 7T°. 



(50) 
(51) 



(45) 



and in the non-relativistic limit (m 3 ^> T,p 3 ), this quan- 
tity tends to 



In this example all particles are bosons and we put 
mi = mi for simplicity, which is not quite exact for re- 
action (l5Tj) . In integral (|42j) we have E\ = E2 = m p /2 in 
the p rest frame. The integrand in $>{p) is a symmetric 
function. Then we can write 



®(P3 -> 0) = 2- 



^ — 1 ^ — 1 ^ 

(Tr^^i'ij^-Wii)' (46) $( P p) = 2^ dc —, Jj_ 2 bQ l T - lei ! +bC) — ■ 



(52) 



where 



Decay and production rates in the Boltzmann 
limit 



A useful check of the more complex quantum decay 
case is the Boltzmann limit. We can then omit unity in 
the distribution Eq. (IT5)) . This is possible when 



J m 2 P + P 2 P J 1 - -iml/m 2 p p p 

a= 2T ; b = 2T • (53) 

The integral, Eq. (|52|). can be evaluated in this case as 
2T 2 



» 1, 



(47) 



HPp) = 



b(T^ 2 e 2a - 1) 



that is, when Tj <C 1 or T <C m^/2. The condition 
T <C m 3 /2 comes from the fact that the minimal energy 
of lighter particles is m%/2 in the particle 3 rest frame. 
In this limit the decay time in the particle 3 rest frame 
from Ea.(|lD|) t' -> r so that from Eq.([2S} we have, for 
the average decay rate r in the reference frame (the rest 
frame of the bath) , 



x fc + ln 1+ M .(54) 



To 



f °°p 2 dpe E ^ T 
j °° p 2 dpe E ^ T m 3 /E 3 

K 2 ( mi /T) 



(48) 



(49) 



Equation (|3"9"]l shows that the average decay time r 3 in 
the laboratory frame is proportional to the (inverse) aver- 
age of the Lorentz factor of particle C3 . We address this 
effect next in a quantitative manner; the ratio of T3 to to 
is shown in Fig|4] as the dotted line. For T <C m 3 this 
ratio goes to unity because the Lorentz factor becomes 1. 
For large T, the rate increases because of the higher aver- 
age energy of particle C3, that is, the increasing average 
Lorentz factor 7. Therefore, for the low density classical 
limit with Tj 1 the average particle life time increases 
with T owing to relativistic effects. However a different 
result can arise for a dense quantum medium. 



Then we substitute $ into Eq. (|3"9"| and using Ea.([2"5]l we 
can calculate p decay and production rates. To calculate 
T3 — > t we use definition (f2U|) . 

In Fig. |4] we present the p decay time in the labo- 
ratory frame normalized by its decay time in the rest 
frame in a vacuum as a function of temperature T for 
T p = T T = 1, [solid (red) line], T p = T^ = 2.56 (dashed 
line), T p = T n = 0.1 (dash-dot line); the dotted line 
shows the Boltzmann limit, Eq. (|4"9")l . We consider the 
range of temperatures between 50 and 150 MeV, which 
includes the QGP hadronization temperature (~ 140 - 
165 MeV). 



We show the case T, 



0.1 (dot-dashed line) to 



check the transition to the Boltzmann limit. We can see 
that for this case, the result is close to the Boltzmann 
approximation for our range of T, as expected. In the 
case of chemical equilibrium, T p = T K = 1, the solid line 
in Fig.@]shows a relatively small, 10-15% increase in lifie 
span. We finally consider a supersaturated pion state, 
T p = Ttt = 2.56, which can arise after supercooled QGP 
hadronization near T = 140 MeV [||. For small T < 
mp/2, the ratio t/tq is near the Boltzmann limit, close to 
unity, because for such a small T, when the Boltzmann 
limit is applied, decay time r does not depend on T. 
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FIG. 4: (Color online) The ratio t/to as a function of tem- 
perature T in the reaction p -H- nn. The [dotted (blue) line] 
is for the Boltzmann limit showing only time dilation. Near 
this limit [dash-dotted (green) line] dilute system T p = T n — 
0.1. Solid (red) line and dashed (turquoise) lines represent 
T p = T,r = 1 and T p = T n = 2.56, respectively. 

When T increases quantum the effects dominate and r 
decreases with increasing T . In general, the larger T is, 
the more rapidly r decreases with temperature. 

Here we do not consider in depth the p-meson density 
evolution in heavy-ions collisions, because without doubt, 
our limited system (just a few hadron states) is not suf- 
ficiently realistic to capture the physics of the p in dense 
matter. Moreover, the p yield is a probe of the hadron 
density temporal evolution and thus still more difficult 
to describe precisely. This can be seen as follows: In the 
Boltzmann low-density limit, the chemical equilibration 
time is T p 1.7 fm. Our result shows that the pion high 
density quantum medium effects causes an increase in p 
width, that is a decrease in equilibration time to r p s» 1.25 
fm, accelerating p-meson chemical equilibration near the 
hadronization temperature. The kinetic-phase time scale 
in heavy-ion collisions, when hadrons interact is near 2-3 
fm [||. This means that the p-meson chemical evolution 
is dependent on the ambient hadron density and thus is 
intricately connected with the dynamics of fireball ex- 
pansion. 

2. <f) meson evolution in heavy ion collisions 



tion, individual hadrons can continue to rescatter into 
resonances in what we call the kinetic evolution phase 
or thermal hadronic gas. This scattering effect does not 
materially change the final stable particle yields, but it 
affects the yields of resonances observed by the invariant 
mass method. The temperature of kinetic-phase freeze- 
out is expected to be near 100 MeV. After kinetic freeze- 
out, hadrons expand without interactions, via decay only. 

The 0-meson has, on the hadron reaction scale a rel- 
atively small width ss 4.26 MeV, that is rj 46 
fm, which is much longer than the duration of the ki- 
netic phase. About 83% of 0-mesons decay into K + K. 
Therefore we consider here <f> evolution in the reaction: 

K + K o <t>. (55) 

We do not consider here the decay channel 4> — » p + 7r, 
which is about 15% and can influence our result at this 
level. Moreover, the <j> inelastic scattering in two-to-two 
particle reactions also has a noticeable influence on <f> 
yield, about 15% suppresion [Io|. In [Io| only the (j) decay 
was included, without the reverse reaction, assuming an 
initial equilibrium yield at hadronization. We show here 
how the reverse reaction and nonequilibrium hadroniza- 
tion conditions can influence the resulting <fi yield. The 
effect from the full one-to-two reactions, Eq. (|55"|) . can be 
added to that from two-to-two particles reactions. 

We did a similar study previously for baryon reso- 
nances A(1232) and £(1385) [? ]. We found that in 
the case of initial non-equilibrium yield, when we have 
an overabundance of stable particles at hadronization, 
the resonance production can be greater than the reso- 
nance decay. Decay becomes dominant when the temper- 
ature drops with expansion, and the lighter mass decay 
product states becomes statistically preferable. The fi- 
nal resonance yield depends on the study of the balance 
between these two effects. 

The </>-meson width is smaller than the width of 
these baryon resonances, and its yield change during 
the posthadronization kinetic phase is expected to be 
smaller. However, the rather low threshold energy, 
— rriK — tuk ~ 30 MeV, could mean that cf> produc- 
tion is dominant over a longer period of time than in the 
aforementioned case of baryon resonance. The purpose 
of this short study is to determine how much the yield 
of (j> can change during the kinetic phase owing to kaons 
fusion compared to its yield at hadronization. We do not 
study how relativistic and quantum effects influence the 
reaction, Eq. (1551) . relaxation time, because for the range 
of temperature considered these effects are small. 

Considering that the mass of all particles involved is 
greater than the temperature it is possible to use the 
Boltzmann distribution for <fi and K: 



We consider here (/>-meson yield evolution in a thermal 
hadronic gas after QGP hadronization formed in heavy- 
ion collisions. The temperature of QGP hadronization 
can be within the range 140-180 MeV. After hadroniza- 



T 3 

2^3 



V 

Nk • < j. 



(56) 
(57) 
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FIG. 5: (Color online) The yield ratio <j>/<j> for hadronization 
temperatures T = 140 MeV [solid (blue) line], T = 160 MeV 
[dashed (green) line] and T = 180 MeV [dash-dotted (red) 
line] as functions of ambient temperature T. 



where Xi = nii/T and 1^2 0*0 is the Bessel function. We 
proceed as in Ref. @, using Eq. (|3Tj) . 

Initial conditions in the kinetic phase are defined by 
conditions at QGP hadronization. We introduce the 
initial hadron yields in a framework of a rapid QGP 
hadronization with all hadrons produced with yields gov- 
erned by the entropy and strangeness content of QGP by 
quark recombination. In this model the yields of mesons 
and baryons are controlled by the constituent quark fu- 
gacity 7 g : 



-y-0 



T 

J- A. 



4> = ls 



(58) 



Thus for 7 g > 1 we have the condition < Tk^k- 
At first, the reaction goes toward <f> production until the 
4> density reaches the equilibrium point when the right 
hand side of Eq. (pTTj) is 0. If the 4> density has enough time 
to reach this point, it begins to decrease again because 
the temperature decreases owing to expansion. 

For each entropy content of the QGP fireball, the cor- 
responding fixed background value of j q can be found 
once the hadronization temperature is known Q. For 
T = 140 MeV pions form a nearly fully degenerate Bose 
gas with 7 g ~ 1.6. In the following discussion, besides 
this initial condition, we also consider the value pairs 
T = 160 MeV, 7 g = 1.27 and T = 180 MeV with j q = 1. 
The value of 7 S > 1 plays no significant role as, in the re- 
action considered, Eq. (|55[) the number of strange quarks 
present is the same. 

In Fig. [S]we present results for the ratios 4>/<fio at differ- 
ent hadronization temperatures as functions of temper- 
ature T, beginning from the presumed initial hadroniza- 



tion temperature T through T min = 90 MeV. (jP is the 
initial yields obtained at each hadronization tempera- 
ture. For hadronization temperatures To < 180 MeV 
(j q = 1.6), we initially have < Tr-Tr-. In these cases 
the master equation leads to an initial increase in the 
yield of resonances. In the case To = 140 MeV, when the 
effect is largest, this increase in (f> yield continues over 
the full range of temperature considered. However, the 
effect is relatively small, about 7%, owing to the small 
<f> width. For hadronization temperature T — 160 MeV, 
when 7 9 = 1.27 is smaller, the increase in yield is smaller, 
and at T = 105 MeV the cj> yield begins to decrease slowly 
owing to the dynamics of the expansion. 

We note that for T > 180 MeV there is always a slow 
depletion of the cf> resonance yield. This result implies 
that the observed yield of <j> has a systematic +7%/-4% 
uncertainty due to kaon rescattering in the medium. For 
comparison in flo| the effect from cj) decay only for equi- 
librium yield at hadronization (j q = 7 S = 1) was de- 
termined to be -7.5% and the effect from 2-to-2 parti- 
cles reactions was -15%. Therefore <f> production in kaon 
fusion for noncquilibrium hadronization conditions may 
have an enhancement effect of about 15% on the final 
<f> yield, compared to the scenario where the 4> can only 
decay after in-cquilibrium hadronization formation. 

We cannot compare with experimental results consid- 
ering only the kaon fusion reaction, as it was argued in 
Ref. [10( that certain two-to-two reactions and possibly 
other processes can influence the yield. We note only 
that kaon fusion can add to the observed (j> yield [Til ]. 



B. Freeze-out processes in the early Universe 

1. 7T° at T < 

As mentioned in Sec. I, it is interesting to examine the 
mean life time of ir° in the end of the hadronic-gas stage 
of the Universe where the temperature drops to a mega- 
electron- volt level in the low teens. Then the reaction 

7T° <-> 7 + 7 

determines the abundance of tt°. 

The difference from the previous example is that the 
photons are massless and they are always in chemical 
equilibrium in the early Universe (Ti = T2 = 1). Then 
we can rewrite function (|45p as 

e (b-a) _ g-(a+6)> 

$Otto) = , , 0." ~ ( b + In ( 1 + ' 



b(e 2a - 1) 



with 



irr 



-Plo 



2T 



1 ~e b 



6 = ^. 



2T 



(59) 



(60) 



Again we use Eq. (|59"|) and (j2"5j) we can calculate n decay 
and production rates. To calculate r for 7r° decay in 
matter we use definition (|2"9")l . 
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FIG. 6: (Color online) The ratio t/t~o for n° decay /production 
as a function of temperature T. Dashed (blue) line is for 
a dilute system Y^o = T 7 = 0.01 (Boltzmann limit); solid 
(green) line is for a thermal chemically equilibrated system 
Y^o = 1. 

In figure [5] we show ratio of 7r° decay time T3 — > r in 
the presence of thermal particles to the decay time in 
vacuum in ttq rest frame: t„o/t^ - In this figure a wider 
range of temperature is shown 1 — 200 MeV. For T^o = 1 
the ratio r^o /r° the temperature dependence is similar 
to that for p decay, considered in previous chapter. It 
increases at first due to relativistic time dilution effects. 
Then, after T w 20 MeV, r goes down slowly with tem- 
perature, when the quantum in-medium effect becomes 
important. Quantum in-medium effects arise here mostly 
from photons. They compensate relativistic Lorentz fac- 
tor effect when T is about m^; compare the lines in Fig. [5] 
for = T 7 = 0.01 [dashed (blue)] and T T = T 7 = 1 
[solid (green)]. Note that when the yield of pions is small, 
that is, only is small, the result is almost the same as 
in chemical equilibrium, = T 7 = 1. 

As long as the ttq reaction relaxation time r is much 
shorter compared to the Hubble expansion time T/T 
1 If. there is chemical equilibrium in the Universe with 
T^n = 1. Freeze-out from chemical equilibrium arises 
when condition Eq. (j3"6")l is satisfied. Because r rj to = 
8.4 10 _17 s the condition Eq. (|36|) is always satisfied where 
ir° can exist. Only at unrealistically large temperatures 
this condition can be violated. 

Therefore we conclude that for the temperature range 
of interest, between few MeV and 180 MeV, the ir° are in 
chemical equilibrium with photons because of their fast 
reaction rate. Note that weak interaction process such as 
neutron decay n — > p + e~ + u e is 20 orders of magnitude 
slower, and the Universe expansion rate can dominate the 



neutron decay rate e.g. at T > 0.1 MeV, before having 
a good chance to decay, neutrons are thus available to 
enter nuclear reactions. 

The pion, and muons equilibrium density is large un- 
til temperatures near a few MeV, and as a result, they 
participate in reactions with each other, nucleons and 
other particles; for example see the next sections. The 
importance of this realization of pion chemical equilib- 
rium in the Universe is that it implies that all hadron 
species, driven by pions also remain in chemical equilib- 
rium. Their abundances can thus be computed using the 
chemical equilibrium hypothesis, as done in Ref.[4]. 

2. tt , fj,^ , and u, v equilibration /freeze-out 

In the laboratory the dominant tt^ production reaction 
is the pion charge exchange reaction, Eq. ([S]), which we 
considered in Ref. Q . These reactions also can take place 
in the early Universe. However, given the much slower 
evolution of the early Universe, we also now encounter 
reactions involving neutrinos, Eq. ©; the related in- 
vacuum, weak-decay lifespan of the 7r ± is to = 2.60 x 10~ 8 
s. 

In Fig. [7] we show the relaxation times in units of tq for 
equilibration, Eq. (I2U1) . in reaction ©, as functions of 
temperature: near T ~ 160 MeV the lifespan is enhanced 
by a factor of 3 for thermal equilibrium densities with Ts 
=1 [solid (blue) line in figure [7] mostly owing to Pauli 
blocking of the decay products. The time dilation owing 
to thermal motion, which also prolongs the life span, has 
a smaller effect, visible in the Boltzmann limit, which we 
study for a dilute system with Ts=0.01, [dashed (green) 
line] . 

Interestingly, as we next show, the process, Eq.©, is 
the fastest mechanism of neutrino equilibration in a wide 
range of temperatures relevant here, T > 7 MeV, but the 
f^-freeze-out condition is at a lower T and seems to be 
controlled by the reaction [ljj 

e + + e~ O v e ^ + De,^ (61) 

which we also consider now. The neutrino oscillation 
effect assures that all neutrinos remain in equilibrium as 
long as one is strongly coupled to the system. 

In Fig. [8] we show muon-neutrino equilibration time 
in reaction Eq.© (solid line, blue); recall that to ob- 
tain this relaxation time from the results shown in Fig. 
[7J we need to replace the ir density in the nominator of 
Eq. (l29]) by the density of v. This relaxation time in- 
tersects the Universe expansion rate at T ps 5.5 MeV. 
Freese et al [l2[ obtain the relaxation time as a function 
of T assuming the neutrino chemical potential /i„ <C T 
in reaction Eq. ([BTj) : 

7i,„ (ee) - (0.1G F T 5 )-\ r Mee) = (0.6G F T 5 )-\ (62) 

where Gf = 1.1664 10 -5 GeV -2 is the Fermi constant. 
These two results are shown in Fig. [5] We see that the 
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FIG. 7: (Color online) ir ± relaxation time as a functions of T 
in reaction Eq.©, in thermal equilibrium [solid (blue) line] 
and in the Boltzmann limit, obtained for T =0.01 [dashed 
(green) line]. 



FIG. 9: (Color online) Equilibration times as functions of 
temperature, for ty° (solid line, red), tt [dash-dot (green) 
line] , jti* (dashed line, blue) and tt ~ T/(Hm 7T ) [dotted 
(turquoise) line]. 




T [MeV] 

FIG. 8: (Color online) Relaxation time for neutrino equili- 
bration as function of temperature compared to the Universe 
expansion time 1/H [dashed (turquoise) line]. Solid (blue) 
line for reaction Eq.(]6} with equilibrium densities (T =1); 
dash-dotted (green) line , and dotted (red) line are for reac- 
tion (|61|l for muon and electron neutrino, respectively. 



muon-neutrino freeze-out temperature according to reac- 
tion ([6]) is slightly higher than that according to reaction 
(|61|) . The temperature of the neutrino decoupling in re- 
action Eq. dH]) is T d =! 3.5 MeV for and T d = 2.0 MeV 
for v e . 

For a wide range of temperatures, to as low as 7 MeV 
neutrino chemical equilibration by reaction ([6]) is dom- 
inant. This example shows that reactions with chemi- 
cally equilibrated pions and muons can have an influence 
on other, even lighter particle evolution for temperatures 
T«m. 

Muons can be equilibrated by reaction ([7]) and by the 
1 <-> 3 reaction [Eq.®]. We do not consider this type 
of reaction in detail here. For low temperatures, T <C 
nifj,, when relativistic and medium effects are small, we 
assume that the muon decay time and reverse reaction 
relaxation time are nearly the muon lifespan in vacuum, 
T = 2.201Cr 6 s. 

In Fig. [S] we show relaxation times for dominant re- 
actions for pion and muon equilibration. For re- 
action, Eq.©, becomes dominant over reaction ([5]) at 
T « 6 MeV. For the reaction © becomes dominant 
at T sa 4 MeV. Therefore at these low temperatures, re- 
laxation times for ^ and 7r equilibration becomes con- 
stant and far below the Universe expansion rate and tt 
[dotted (turquoise) line]. We conclude that ^ and 7r 
stay in chemical equilibrium. This does not mean that 
they play an important role in the global physics of the 
early Universe, because just at these temperatures muon 
and pion densities begin to drop rapidly and soon their 
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yield is negligibly small, far below the nucleon density in 
the Universe. 



V. CONCLUSIONS 

We have presented detail of the kinetic master equation 
for a process involving the formation of an unstable par- 
ticle through reaction ([TT]) in a relativistically covariant 
fashion. Assuming that all particles in the process are in 
thermal equilibrium, we calculated the thermal averaged 
decay and formation rates of the unstable particle. Us- 
ing the time reversal symmetry of quantum processes, we 
have shown that the time evolution of the density of an 
unstable particle is given by (fT5|) . Therefore in chemical 
equilibrium the particle fugacities are connected by (|18p . 
We have explicitly derived the thermal decay rate of an 
unstable particle, obtaining Eq. (|44l) . 

The general properties of the thermal particle de- 
cay/production kinetics have led us to consider the relax- 
ation time defined by Eq. (j2"9")l . which results in a greatly 
simplified kinetic equation Eq. (|31[) . The medium modifi- 
cation of reaction rates we encountered are all caused by 
final-state quantum effects, Bose enhancement, and/or 
Fermi blocking, absent in the Boltzmann limit. More- 
over, we note the presence of kinematic effects, in that 
all lifespans of particles are time dilated owing to their 
motion with respect to the thermal bath rest frame. 

In the present formalism, we assumed that the decay 
width of an unstable particle is much smaller than the 
temperature T. This approximation is safe in the exam- 
ples we have discussed above, except perhaps the case 
of p decay, where some corrections may be needed. For 
the formation of heavy resonances, whose decay width 
becomes appreciable compared to the temperature T, we 
may need to include the finite-width effect on the mass 
of an unstable particle in the thermal distribution. Such 
effects on the statistical partition function and the equa- 
tion of state of a system have been studied based on the 
virial expansion method [l3j . The correction for a kinetic 
equation in such a case has also been studied fl4j . 

We have presented several examples, p <h- tt + n, 
cj> K + K, 7T° o 7 + 7, and 7r ± O p ± + i^(tv), and 
explored the physics cases of hot hadron matter created 
in laboratory heavy ion reactions, and the early Universe 
from the condition of hadronization down to the temper- 
ature of several mega-electron volts. The two first pro- 
cesses can take place in both circumstances. The third 
process is important to the understanding of how the 
hadronic fraction evolves with the expansion of the Uni- 
verse. The last process we considered appears to be the 
dominant mechanism of neutrino equilibration over the 
entire temperature range, except close to neutrino freeze- 
out, a result that requires further refinement allowing for 
finite chemical potentials. This example also shows that 



heavy (m >> T) chemically equilibrated particles can 
be important in the evolution of other particles includ- 
ing lighter more dense particles yields at relatively low 
temperatures. 

The equilibration-relaxation time for 7r° decay remains 
close (within 25%) to the relaxation time in vacuum for 
a large temperature range. This occurs because the rel- 
ativistic effect (Lorentz factor) is compensated by the 
quantum medium effect. This time is short compared to 
the Universe expansion time for all temperatures of in- 
terest here, below the QGP hadronization temperature, 
when 7T° hadrons are created. Therefore ir° always stays 
in chemical equilibrium with radiation for the tempera- 
ture range of interest. 

As long as 7r° is abundant it can participate in reac- 
tions with other hadrons and influence the dynamics of 
the Universe evolution. Here we also considered tt^ evo- 
lution, in the fourth reaction given and their interaction 
with 7T°. We showed that pions and muons (mesons) stay 
in chemical equilibrium throughout the evolution of Uni- 
verse, despite their large mass. They can be involved in 
reactions with nucleons, a topic we postpone to a future 
study, down to temperatures where the meson density 
drops well below the nucleon density. The contribution 
of mesons disappears from the entropy and the degener- 
acy g only at the relatively low T « 10 MeV (see Fig. [3]). 

Our study of the <f> evolution in thermal hadron 
medium after QGP hadronization in heavy ions colli- 
sions, suggests a possible slight modification of the ob- 
served <j) yield: compared to initial production, an in- 
crease in hadronization at T — 140 MeV, 7 g = 1.6, at 
the level of about 6%-7%, or a suppression of about 4% 
for hadronization at T — 180 MeV, j q = 1. 

To conclude, we have presented here the process of 
decay and and re-creation of unstable particles, and 
studied special cases of relevance to heavy-ion collisions 
and the early Universe. Our results indicate that the 
early Universe was in chemical equilibrium throughout 
its evolution and that the first freeze-out occurs when 
neutrinos decouple. 
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